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The recent discovery of superconductivity (SC) in MgB2 |lj] stimulated a significant interest in this system. One 
of the first questions is whether MgB2 represents a new class of superconductors, or whether it may be related to 
other known classes in terms of its bonding and electronic properties. Is it a unique system or just a representative 

^vj ' of some family of compounds with similar or even better properties? The latter case seems unlikely after all 
studies [|[ H of possible dopings such as Mgi_2;Ala;B2, MgB2_a;C:E, Mgi_j;Lij;B2 and Mgi_2;Mn:rB2 showed loss of 

C^^ ' SC, and many attempts to dope MgB2 with other elements failed. 

^— ^ [ The crystal structure of MgB2 may be regarded as that of completely intercalated graphite H with carbon 

replaced by boron, its neighbor in the periodic table. In addition, MgB2 is formally isoelectronic to graphite. 
Therefore, chemical bonding and electronic properties of MgB2 are expected to have strong similarity to those of 
C^ ' graphite and graphite intercalation compounds (GIC's), some of which also exhibit SC. 

The search for high-temperature SC in carbon phases started in 1965 when SC with Tc — 0.55K in KCs was 
reported and subsequently explained in terms of the interaction of tt bands with the longitudinal phonon modes 

' ^ , of the carbon layer H. Similar conclusions for the electronic structure at the Fermi level were derived for LiCe. 
The highest Tc achieved for a GIC was 5 K g. A parallel development was that of SC in Bechgard salts, the 
(J , organic charge-transfer systems. Within this family of materials exhibiting carbon n band conductivity, Tc was 
raised to 12.5 K 0. It should be noted that the majority of these compounds are unstable at normal conditions 
and high-pressure experiments are most common. 

j^ ' Later research shifted to the area of carbides of transition metals where metallicity and conductivity are mostly 

H I due to transition metal atoms. Only in 1991 was it found that alkali-doped fuUerenes also exhibit SC with the 
highest Tc of 33 K g . This group of 3D carbon-based metals has a modified (compared to graphite) coupling of 
lattice vibrations to Fermi electrons due to the curvature of the molecule, but the conducting states still derive 
from the graphite tt band. 

In this paper we will show that in spite of the structural similarity to intercalated graphite and, so some extent, 
to doped fuUerenes, MgB2 has a qualitatively different and rather uncommon structure of the conducting states 
setting it aside from both these groups of superconductors. The peculiar and (so far) unique feature of MgB2 is 
the incomplete filling of the two a bands corresponding to prominently covalent, sp^-hybridized bonding within 
the graphite-like boron layer. The holes at the top of these a bands manifest notably two-dimensional properties 
and are localized within the boron sheets, in contrast with mostly three-dimensional electrons and holes in the 
TT bands which are delocalized over the whole crystal. These 2D covalent and 3D metallic-type states contribute 



almost equally to the total density of states (DOS) at the Fermi level, while the unfilled covalent bands experience 
strong interaction with longitudinal vibrations of the boron layer. 

Below we discuss the properties of MgB2 following from its theoretical treatment in the local density approx- 
imation (LDA). However, an unusually strong non- linear electron-phonon coupling (EPC) M together with the 
probable non-adiabaticity ||l^ raise a question whether the ground state and the excitation spectrum may be 
adequately treated in LDA. In this situation it is especially important to check the LDA predictions against the 
experimental data related to the electronic structure of MgB2 and doped alloys. Therefore we also try to provide 
the basis for such comparison using both our own results and those available in the literature. 

This paper is organized as follows. The electronic structure of bulk and surface of MgB2 and its relation to 
GIC's is described in Section 2 along with the effects of doping. Here we justify our choice of the Mgi_2:Ala;B2 
system for detailed studies allowing one to analyze trends associated with band filling and their relation to the loss 
of SC. We believe that such analysis is more valuable as opposed to the studies of a single system due to many 
uncertainties (physical and numerical) of current band structure calculations and experiments. The relevance of 
the theoretical DOS at the Fermi level N = N{Ep) as well as the anisotropy of conducting states are analyzed in 
Section 3 in connection with NMR data for Mgi_a;Ala;B2 alloy. Short description of the EPC studies is presented 
in Section 4 along with the comparison of theoretically predicted N with that deduced from the experimental data. 
Optical calculations for MgB2 and AIB2 single crystals are described in Section 5. 

2 Electronic structure and bonding 

MgB2 occurs in the AIB2 structure. Boron atoms reside in graphite-like (honeycomb) layers stacked with no 
displacement [Q forming hexagonal prisms with the base translation almost equal to the height, a = 3.085 (3.009) 
A and c/a — 1.142 (1.084) for MgB2 (AIB2). These prisms contain large, nearly spherical pores occupied by Mg 
atoms. As in graphite (i?intra=l-42 A), the intralayer B~B bonds are much shorter than the interlayer distance, 
and hence the B-B bonding is strongly anisotropic. However, the intralayer bonds are only twice as short as 
the interlayer ones compared to the ratio of 2.4 in graphite, allowing for a significant interlayer hopping. For 
comparison, the interatomic distance between nearest neighbors is 1.55 A in diamond and 1.4-1.45 A in the Cgo 
molecule. 

Early studies of the band structures of MgB2 and AIB2 [12-15], as expected, showed their notable similarity 
to that of graphite (see e.g. M and references therein). Full-potential band structure calculations of AIB2 are 



available in the literature |15 . A recent paper |17| reported the results of the studies of MgB2 and AIB2, but the 
structure of conducting states was not addressed. 

Below we discuss the electronic structure of MgB2 and some related compounds obtained [18-20] using the 
Stuttgart TB-LMTO (ASA) code, full-potential LMTO (FLMTO) Q and full-potential LAPW (FLAPW) meth- 
ods. It was found that the addition of gradient corrections to the local exchange-correlation potential has only a 
tiny effect on the results. It also appears that a general analysis of energy bands for MgB2 may be safely done in 
ASA, while a reliable treatment of charge densities and anisotropy of transport properties requires more accurate 
full-potential calculations. The band structure of AIB2 is much more sensitive to the choice of atomic spheres in 
ASA, and full-potential treatment is imperative. Other methods were also used to calculate the electronic structure 
of MgB2 with very similar results. 

The energy bands, DOS and the Fermi surface of MgB2 are shown in Figs. ^ |p and ||a. As expected, the 
bands are quite similar to those of graphite with three bonding a bands corresponding to in-plane spxPy (sp^) 
hybridization in the boron layer and two tt bands (bonding and antibonding) formed by aromatically hybridized 
boron pz orbitals. Both a and tt bands have strong in-plane dispersion due to the large overlap between all p 
orbitals (both in-plane and out-of-plane) for neighboring boron atoms. In particular, the total width of the tt 
bands is 17 eV (which gives ^ ^* — 0.14, see below). The interlayer overlaps are much smaller, especially for p^y 
orbitals, so that the k^ dispersion of a bands does not exceed 1 eV. From the other hand, in contrast to graphite 
and GIC's, two of the a bands are filled incompletely, as it was first noted in Ref. [Q. Together with weak kz 
dispersion this results in the appearance of two nearly cylindrical sheets of the Fermi surface (see Fig. ||a) around 
the F-A line. As we will see below from the analysis of the charge density distribution, these unfilled a bands with 
boron p^y character fully retain their covalent structure. Conducting covalent bonds represent a peculiar feature 
of MgB2 making it an exotic compound probably existing on the brink of structural instability. The blue (red) 
hole-type (electron-type) tubular network in Fig. 0a corresponds to the bonding (antibonding) tt band. 
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Figure 1: Energy bands of MgB2 from Rcf. ||l8[. The radii of black (red) circles are proportional to the boron p^y 
(pz) character. 
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Figure 2: Total density of states for (a) MgB2 and (b) AIB2. Zero energy corresponds to the Fermi level. 







Figure 3: Fermi surfaces of (a) MgB2; (b) AIB2; and (c), (d) MgB2 with E2g distortions of ±0.07ao (see text). 



It is seen in Fig. ||a that the tt bands form two planar honeycomb tubular networks: an electron-type sheet 
centered at A;^ = (red) and another similar, but more compact, hole-type sheet centered at kz = tt/c (blue). These 
two sheets touch at some point on the K-H line. Note that the hole-type sheet is close to the electronic topological 
transition (ETT) at the M point corresponding to the breakdown of the tubular network into separate shell-like 
pockets (at 0.25-0.30 eV above Ep). Although the singularity in the electronic properties is weak, the proximity 
of this ETT results in a strong coupling with the in-plane £^2g phonons, and vibrations of moderate amplitude are 
able to break the links at the M point (see the discussion of Figs. 0c and Hd below). In addition, this ETT takes 
place with electron doping (see below). 

From Fig. 0a it is clear that the Fermi surface of MgB2 has characteristic features (the cylindrical sheets and the 
tubular links at M and L points) that almost completely determine its topology. These features are associated with 
relatively small electronic orbits that are especially suitable for identification in de Haas-van Alphen experiments. 
Such experiments on single crystals are highly desirable for the verification of the band structure obtained in LDA. 

In order to examine the relation between the band structure of MgB2 and that of graphite in more detail we 
studied the following hypothetical sequence of intermediate materials: carbon in the 'primitive graphite' (PG) 
lattice with no displacement between layers as in MgB2, using graphite lattice parameters; boron in the PG lattice 
with a as in MgB2 and c/a as in graphite; boron in the PG lattice with a and c/a as in MgB2; LiB2 in the same 
structure; MgB2 itself. The results of some of these calculations are shown in Fig. |[ 

The band structure of PG carbon shown in Fig. ^ is very similar to that of graphite |Q with the appropriate 
zone- folding for a smaller unit cell. (This is quite natural because of the weak interlayer interaction.) Boron in 
the same lattice dilated to match the MgB2 in-plane lattice parameter (not shown) has nearly identical bands 
with the energies scaled by the inverse square of the lattice parameter, in agreement with common tight-binding 
considerations ||2^ . Fig. |jb shows the natural enhancement of the out-of-plane dispersion of the tt bands when 
the interlayer distance is reduced. Figs. ^ and ^ demonstrate that 'intercalation' of boron by Li or Mg produces 
a significant distortion of the band structure, so that the role of the intercalant is not simply one of donating 
electrons to boron's bands (which would return the band structure to that of PG carbon shown in Fig. ^). The 
main change upon intercalation is the downward shift of the n bands compared to a bands. For Li this shift of 
~1.5 eV is almost uniform throughout the Brillouin zone. Replacement of Li by Mg shifts the tt bands further, but 
this shift is strongly asymmetric increasing from ~0.6 eV at the F point to ~ 2.6 eV at the A point. In addition, 
the out-of-plane dispersion of the a bands is also significantly enhanced. In LiB2 the filling of the bonding p,j:y 
bands is nearly the same as in PG boron, while in MgB2 the Fermi level shifts closer to the top of these bands. 

The lowering of the tt bands in MgB2 compared to PG boron is due to stronger interaction of boron p^ orbitals 
with ionized magnesium sublattice compared to Pxy orbitals. This lowering is greater at the AHL plane compared 
to the FKM plane, because the antisymmetric (with kz — tt/c) overlap of the boron's p^ tails increases the electronic 
density close to the magnesium plane where its attractive potential is the strongest. 

The nature of bonding in MgB2 may be understood from the charge density (CD) plots shown in Fig. pi As it is 
seen in Fig. ^a, bonding in the boron layer is typically covalent. The CD of the boron atom is strongly aspherical. 



and the directional bonds with high CD are clearly seen (see also Ref. 1 17 ) . The CD distribution in the boron 
layer is very similar to that in the carbon layer of graphite |16[| . This directional in-plane bonding is also obvious 
from Fig. ^d showing the CD in the cross section containing both Mg and B atoms. However, Fig. |^b also shows 
that a large amount of valence charge does not participate in any covalent bonding, but is rather distributed more 
or less homogeneously over the whole crystal. Further, Fig. |^ shows the difference of the CD of MgB2 and that 
of hypothetical NaB2 in exactly the same lattice. Not only does it show that one extra valence electron is not 
absorbed by boron atoms but is rather delocalized in the interstitials; it also shows that some charge moves outward 
from boron atoms and covalent in-plane B-B bonds. Fig. |^d shows the CD difference between the isoelectronic 
compounds MgB2 and PG carbon (C2). In MgB2, the electrons see approximately the same external potential as 
in C2, except that one proton is pulled from each C nucleus and put at the Mg site. It is evident that the change 
C2— >MgB2 weakens the two-center a bonds (the charge between the atoms is depleted) and redistributes it into a 
delocalized, metallic density. 

A numerical reconstruction of the electronic charge density from the synchrotron radiation data for a powder 
MgB2 sample [p4J supports this general picture. The charge density found for 15 K is, in fact, very similar to 
our Fig. ob and shows all the important features discussed above including the distinct covalent bonds within the 
boron sheets, the strongly ionized Mg, and the delocalized charges in the interstitials. Further, the Fourier maps 
obtained pl| for the single crystals also clearly show the covalent sp^ hybrids in the boron layer and no covalent 
bonding between B and Mg atoms. 

Because of the coexistence of two different types of conducting states, it is desirable to find the contributions 
to the total DOS and transport properties from separate sheets of the Fermi surface originating from 2D covalent 
and 3D metallic-type bonding. Such decomposition is shown in Fig. H for the total DOS (see also Ref. P5|) and 




Figure 4: Band structures of: (a) top left: primitive (AA stacking) graphite (PG), a = 2.456A, c/a — 1.363; (b) 
top right: PG boron, a — 3.085A, c/a = 1.142 (as in MgB2); (c) bottom left: LiB2 in MgB2 structure, same a and 
c/a; (d) bottom right: MgB2, same a and c/a. Energy is in eV relative to Ep. The order of occupied bands in the 
F point is a bonding with boron s character, n bonding with boron pz character, and a bonding with boron p^y 
character (double degenerate). 




Figure 5: Pseudocharge density contours obtained in FLMTO. The unit cell is everywhere that of MgB2. Darkness 
of lines increases with density, (a) MgB2 in (0002) plane passing through B nuclei; (b) MgB2 in (1000) plane 
passing through Mg nuclei at each corner of the figure. B nuclei occupy positions (1/3,1/2) and (2/3,1/2) in the 
plane of the figure. The integrated charge of the unit cell is 8. (c) (1000) plane, difference in smoothed density, 
MgB2 minus NaB2. The integrated charge of the unit cell is 1. (d) (1000) plane, difference in smoothed density, 
MgB2 minus PG carbon. The integrated charge of the unit cell is 0. In (c) and (d), dotted lines show negative 
values. 



for the in-plane (xx) and out-of-plane (zz) components of the tensor aafs = J VaVf3S{s{k) — EF)d'k, where Va is the 
a-component of the electronic velocity. The 3D (metallic-type bonding) and cylindrical (covalent bonding) parts 
of the Fermi surface contribute, respectively, about 55% and 45% to N. N{E) for the hole-type zones rapidly 
decreases with increasing e and already a,t E — Ep « 0.8 eV the total DOS is almost completely determined by 
the 3D electron-type band. The latter contribution is almost constant and probably not related to the change of 
SC properties under pressure or with doping. The corresponding contribution to conductivity exceeds all other 
contributions (more than 50% for azz and Uxx) and is virtually isotropic. Holes in the two cylindrical sheets, as 
expected, have clearly anisotropic behavior contributing nearly 30% to axx and virtually nothing to azz- The 
3D hole-type part of the Fermi surface is also notably anisotropic with predominantly z-axis conductivity. The 
total a has a rather small anisotropy at Ep with Uxx/fTzz — 1.22. The average projections of the Fermi velocities 
v„ = (i;2)i/2 are: Vx = 5.36 • 10^ cm/s, Vz = 4.85 • 10^ cm/s (see also Refs. |l|, |^, |§). The average Fermi 
velocities Vi, — {v'^)J for each Fermi surface sheet v are as follows: internal cylinder, 8.85-10^ cm/s; external 
cylinder, 6.02-10'' cm/s; 3D electronic sheet, 1.12-10^ cm/s; 3D hole sheet, 8.74-10^ cm/s. 

Thus, according to Ref. ||l9|, the structure of MgB2 is held together by strongly covalent bonding within boron 
layers and by delocalized, ^metallic-type^ bonding between these sheets. A peculiar feature of this compound is 
that electrons participating in both of these bond types provide comparative contributions to N. This distinguishes 
MgB2 from closely related GIC's where covalent bonds in the carbon layers are always completely filled, while the 
nearly cylindrical parts of the Fermi surface commonly found in those compounds are formed by carbon-derived ir 
bands which are also responsible for conductivity in pristine graphite |5[| . 

Thin films are of great technological interest; indeed probably the most important application of MgB2 at 
present is that for digital logic in very high speed (> 100 GHz) switching applications. Since the device speed 
scales as the bandgap (and hence with T^,) , MgB2 shows great potential promise as the superconductor of choice. 
These devices are grown as thin films; thus the role of the surface effects is of interest. 

We considered a 7-unit cell (21 atom) slab of MgB2, with 7 cells stacked along the c axis, followed by a larger 
vacuum layer to separate the two faces. A unit cell of MgB2 consists of alternating Mg and B planes along the z 
axis. Thus this slab has two kinds of surfaces: a Mg-terminated surface at z = and a B-terminated surface at 
z = 7c for an ideal, unrelaxed geometry. The slab was fully relaxed; however, it retains the hexagonal symmetry 
of the bulk lattice, so the only allowed relaxations are in the z coordinate (no symmetry-lowering reconstruction 
was considered). A substantial relaxation was found at both surfaces, each contracting inward so as to reduce 
the spacing between layers. The relaxation was most pronounced at the Mg-terminated surface (the surface Mg-B 
spacing contracted by 12%). Fig. |^ shows total DOS (top panel) and partial DOS in sites centered at boron 
spheres, resolved by layer (panels 2-8). The DOS shown was computed for 8000 points in the full Brillouin zone, 
using the electronic structure method of Rcf . [pll . The dotted lines correspond to that of bulk MgB2 : 7 times the 
total DOS of one unit cell of MgB2 (top panel) and B partial DOS (central panel). The first layer corresponds 
to the Mg-terminated surface; the last to the B-terminated surface. Several points are evident. (1) The B partial 
DOS of the central layer is quite similar to that of the bulk DOS, showing that surface perturbations heal in just 
a few monolayers, as expected for a metal. (2) There is a shift in the B partial DOS towards higher energies at 
the B-terminated surface (B7), and to lower energies at the Mg terminated surface (Bl). To a large extent this 
effect is electrostatic: the electrostatic potential at the MT boundary for Bl is about 32 mRy higher than for the 
central layer B4, while the corresponding potential for B7 is about 9 mRy lower. (3) There is an enhancement of 
DOS near the Fermi level at the B-terminated surface. This is in part due to the electrostatic shift, but also there 
is an additional structure in the local DOS corresponding to the formation of the surface band which localizes 
approximately at B6 and B7. (4) There is a slight increase in the total DOS relative to the bulk (top panel), which 
originates in the surface states from the B-terminated surface. 

Now let us discuss the relation of MgB2 to other compounds. The closest existing material is the isostructural 
AIB2 corresponding to the addition of one electron to MgB2. Al is the only neighbor of Mg in the periodic table 
that may be used as a dopant producing an isostructural solid solution with a reasonably wide single-phase region. 
As it is always more useful to study a sequence of similar compounds or doping trends in alloys instead of analyzing 
a single system, below we discuss the trends in the band structure and related properties of the Mgi_2:Ala;B2 alloy. 

Our FLMTO band structure and DOS of AIB2 shown in Figs. ^ and ||b are in excellent agreement with earlier 
results |l|, |l|. 

The Fermi surface for AIB2 is shown in Fig. 0b. The red electron-type network (blue pocket at K) corresponds 
to the antibonding (bonding) n band, while red and yellow pockets at F and H correspond to the 'interstitial' 
electron-type band. In contrast to MgB2, the a bands in AIB2 are completely filled, and the Fermi surface has 
no cylindrical sheets. This difference stems from an extra valence electron, but the rigid band picture does not 



fully describe the effect of replacement of Mg by Al 1 20 . The n bands are pushed further down in the AHL plane 




Figure 6: Rigid band results for MgB2: (a) Total DOS in cV^^/cell; (b) axx and (c) azz in Ry-CQ (for definition of 
(Ta/3 see text), witlr contributions from different sheets of the Fermi surface. Thick solid lines: total; thin dashed: 
3D electronic sheet; dotted: 3D hole sheet; thick dashed: internal cylinder; thin solid: external cylinder. Energy is 
in eV relative to Ep- 
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Figure 7: Total DOS (top panel) and partial DOS in sites centered at boron spheres, resolved by layer (panels 2-8) 
for 7-unit cell slab of MgB2. The dotted lines show total (top panel) and partial boron DOS (central panel) in bulk 
MgB2. The first layer (Bl) corresponds to the Mg-tcrminatcd surface; the last (B7), to the B-tcrminated surface 
(see text). 
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Figure 8: Energy bands for AIB2 



(compared to the a bands), together with the interstitial band |g^ which is unfiUed in MgB2 (~2 eV at the F 
point). In AIB2 the latter interstitial band falls below the Fermi level and forms electronic pockets at F and H 
points. The hole-type tubular sheet of the Fermi surface found in MgB2 shrinks to very small hole pockets at the 
K point in AIB2. 

It has been shown [gj that Mgi_a;Ala;B2 in the < x < 0.4 region forms single-phase solutions at < x < 0.1 
and at 0.25 < x with the same AIB2 structure, while at 0.1 ^ a; < 0.25 the system appears to decompose in two 
isostructural AlB2-type phases with different cja ratios. The Mg-rich phase remains superconducting up to a; ~ 0.1 
where T^ is only 2 K lower than in pure MgB2, while the transition broadens in the two-phase region and there is 
no supeconductivity in the Al-rich phase with x > 0.25 (as well as in pure AIB2 pTl). 

Unfortunately, very little is known about the subtle structural transition in the Mgi_2;Al2:B2 alloy. However, 
as this transition is accompanied by the loss of SC, more detailed experimental characterization of the structural 
transition as well as of the dependence of the electronic structure on doping are highly desirable. In particular, our 
FLMTO calculations show that in the rigid band model at approximately x = 0.17 the band structure undergoes 
an above-mentioned ETT at the M point where the hole-type tubular sheet of the Fermi surface splits into compact 
pockets. In principle, this ETT may be the source of the structural transition, but verification of this connection 
needs further studies. 

Hole doping by replacing Mg by Na or Li seems logical, but the experiments on Mgi_2:Li2:B2 (as well as on 
MgB2_xCa; and Mgi_2;Mn2;B2) revealed the loss of SC Q. Weak SC was observed under pressure in NaC2 |2^ 
which is very unstable. However, it may be possible to form an alloy (Mg,Na)B2 with modest amounts of Na. It 
is evident from Fig. y that the hole-type parts of the Fermi surface change dramatically with band filling. 

The band structure of BeB2 in AIB2 lattice calculated with experimental average or optimized lattice parameters 
is quite similar to that of MgB2 
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|2Cf| . The differences include a notable reduction of the lattice parameters 
a ratio), wider valence bands, larger dispersion of the a bands, and a somewhat different shape 
|. Due to this latter difference, the hole-type and electron-type sheets of the Fermi surface 



(especially of the c 
of the TT bands |g( 

corresponding to bonding and antibonding tt bands become almost identical in BeB2. The ETT at the M point 

(tubular link breaking) shifts approximately 2 eV higher compared to MgB2 and becomes irrelevant for the EPC. 

The hypothetical ZnB2 in MgB2 structure is also very similar to MgB2 in terms of the electronic properties. 

Our calculations for this compound with lattice parameters of MgB2 produced a very similar band structure with 
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Table 1: Partial DOS for s and p orbitals at B site, 10 ^ (eV-spin-atom) ^ 
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Pz 


Px 


MgBa 


3.4 


50 


36 


AIB2 


3.3 


19 


1.9 



a nearly identical Fermi surface. Another recently discovered isoelectronic system LiBC [Q contains both B and 
C. According to our results this system is a perfect insulator, and any substitution of C by B will lead to metallic 
behavior. Experimentally small conductivity was observed in LiBC 129]. Electronic structure of this system is also 
very similar to that of BN |Q . 

3 Nuclear spin-lattice relaxation 

For SC, one of the important parameters is N in the normal state. Several groups have performed LDA calculations 
of this quantity and the results range within 0.68-0.74 states/(eV-f.u.) [18-20, 26]. Experimentally N can be 
determined in many different ways but the nuclear spin- lattice relaxation (NSLR) rate T^ measurements represent 
an excellent opportunity to check experimentally not only the total N and its partial components but also their 
anisotropy, i.e. the distribution between the in-plane and out-of-plane p orbitals. For instance, in the studies of 
N in alkali-doped Cgo the theoretical analysis of NMR data provided a qualitatively new interpretation of NMR 
experiments. 

For NSLR on "B (^ = 2.689/iAr) in MgBz the experimental papers [31-33] reported TTi =180, 155 and 165 
K-sec. The relaxation rates were interpreted in terms of dipolar and orbital contributions due to the low Korringa 
ratio and the known dominance of p states at the Fermi level [18-20, 26]. Below we show the results obtained in 
a more careful theoretical analysis for Mgi„a,Ala;B2 |3J]. As we will see, the studies of NSLR in this system may 
provide valuable information about the genesis of anisotropy of electronic states at the Fermi level. 

We used the following expressions for a monocrystaline material with the hexagonal symmetry ]pq] : 

(Tf 1)^ - {47rkT/h){hjMHFfN^ 
(^r')orb = {^7TkT/h){h^r^H^^^^)' [2Ne'Na'^ - Ne'{Ne' - NA'/)sm' 0] 

for Fermi-contact and orbital contributions, respectively. Here H^^.^^ = 2iXB{r^^)p and Hp = (87r/3)/XB[(ps(0, Ep)^ /4,it 
are the corresponding hyperfine fields, {r~^)p is the expectation value of r~^ over the p state (pj,{r, Ep), and Nr is 



the partial DOS at the Fermi level for the representation F (for more details see Refs. [|35|, |36|). To obtain these 
values we used the LMTO-ASA method. To check the sensitivity of our results to the parameters of calculations 
we used different exchange-correlation potentials and inputs with different radii of the B sphere (both with and 
without empty spheres). The relevant partial contributions to N in MgB2 and AIB2 are listed in Table |^. One can 
see that in MgB2 all p orbitals on the B site have a sizeable DOS, while in AIB2 only pz orbital has a large DOS 
with Npx ~ O.lNpz- The s component in AIB2 becomes relatively more important compared to MgB2 resulting in 
the dominance of the Fermi-contact mechanism of NSLR, as we will show below. In both materials the contribution 
of d states to NSLR is very small. As for the Mg site, the s component of TV is the most important, and we expect 
that NSLR for the Mg nucleus is dominated by the Fermi-contact mechanism. However, below we will focus on 
the ^^B NSLR, because of the boron role in SC. 

To calculate Ti according to a general prescription one has to estimate the values of {r~^)i for different / and 
the electronic density at the nucleus ip'^(0)/4:Tr. We found that the convergence of the total NSLR rate with respect 
to the boron sphere radius in this sp system is worse compared to d metals (from rs = 2.3ao to rs = 2.4ao the 
total T^^ decreases by '^15%). The uncertainty is mostly related to the value of {r~^)pNp. 

In our calculations we used the largest vb that were possible without a significant distortion of the band 
structure, 2.4ao for MgB2 and 2.1ao for AIB2. For these radii we have {{ao/r)^)p ~ 1.11 in MgB2 and 1.37 in 
AIB2. For comparison, the atomic value [B^ for {{ao/r)^)p in B is 0.78. The electronic densities on the nucleus 
aoV's(0)/47r for MgB2 and AIB2 were, respectively, 2.68 and 3.02. 

The contributions to the ^^B relaxation rate for the polycrystalline sample calculated using the general formu- 
las [p3 are given in Table 0. The in-plane and out-of-plane p orbitals in MgB2 have similar DOS, and hence the 
relative magnitude of orbital and dipolar contributions to NSLR is close to the 3/10 rule for p states in a cubic 
crystal described by Obata ^2] • The Fermi-contact contribution is also important and amounts to ^40% of the 
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Table 2: Contributions to (TTi)-i [lO""* (K-sec) 





Contact 


Orbital 


Dipole 


Total 


Enhanced 


Experiment 


MgBa 
AIB2 


12 
21 


30 

1 


9 
1 


51 
23 


81 
26 


56 ]3|, 64 |3|, 61 |3| 



orbital term. The contributions from the d partial waves to the dipole and orbital relaxation rates are small (at the 
order of 1%) due to the low diagonal and off-diagonal densities of states {Nd/NpY ^ 0.02 and {Npd/NpY ^ 0.05. 
The quadrupole contribution to NSLR is negligible due to the rather small ^^B quadrupole moment. 

The values of T^ obtained in such manner correspond to the theoretical 'bare' partial boron TV whereas the 
actual T]~ contains different terms enhanced according to the corresponding susceptibilities. To estimate the 
possible range of enhancement, we calculated the total enhancement coming from the spin susceptibility only 
(which is the case for AIB2). We estimated the effective Stoner exchange parameter / = AE/m — 1.1 eV from the 
splitting of the bands at the T point in the external magnetic field. The corresponding Stoner enhancement of the 
uniform spin susceptibility in the 3D case may be written |Q as 5* = 3/ [(1 — NI) (3 — 2/7V)] w (1 — IN)^"" with 
a K, 1.62, while in the 2D case the enhancement is described by the same formula with a = 2. In our case due 
to the mixed 2D and 3D character of the bands it is not clear what value of a should be used, but the difference 
in the result for a — 1.62 and 2.0 is less than 10%. We used a = 1.9 resulting in the enhancement of T^^ by 
approximately 60% (Table g). The obtained 'bare' and enhanced values provide a range of possible NSLR rates for 
MgB2 which should be compared with the experimental rates [31-33] of (56-64)-10"'* (K-sec)~^. The fact that such 
simple estimate may give a faster relaxation compared to experiments suggests a possible importance of unique 
effects resulting in the lowering of the effective boron N . 

The roles of the three NSLR mechanisms are very different in AIB2 where no experimental data are available. 
According to our theoretical estimation, due to the sharp decrease of the boron p component of N compared to 
MgB2 and its very strong anisotropy (see Table Qj), the orbital and dipolar contributions to NSLR become very 
small, and the NSLR in AIB2 is completely dominated by the Fermi-contact mechanism. This conclusion may be 
verified experimentally by studying the Korringa ratio. The resulting Stoner-enhanced NSLR rate in AIB2 is more 
than three times smaller than in MgB2 (see Table 0) . 

We also calculated the NSLR rate in Mgi_a;Ala;B2 in the rigid band approximation with MgB2 bands. From 
Fig. y one can see how the roles of different mechanisms of NSLR change with doping. The sharp decrease of 
A'' in the 2D sheets of the Fermi surface with doping [|l9| leads to the corresponding lowering of all contributions 
to NSLR, and at the point of the complete filling of these 2D bands we expect a very large Ti. Experimental 
verification of this sharp dependence of Ti on doping in this alloy may be a crucial test of our understanding of the 
electronic structure of this system and is highly desirable. Together with the nuclear quadrupole rcsonanse data 
(which is related to the anisotropy of total charges on different p orbitals) such measurements should help to build 
a general picture of the anisotropy of p orbitals. 

The above calculations have been done for a polycrystalline sample. Because single crystals are becoming 
available, we include our estimations of the anisotropy in the angular dependence of NSLR rate JS^I A + i?sin^ 6. 
For MgB2 we obtained B/A k, —0.06, so that the NSLR is nearly isotropic. The NSLR rate in AIB2 is isotropic 
because it is determined by isotropic Fermi-contact mechanism. 

Summarizing the results of this Section, the orbital NSLR mechanism dominates over the spin-dipolar and 
Fermi-contact ones in MgB2, because the boron p orbitals at the Fermi level are distributed nearly isotropically 
and have a large DOS. Strong anisotropy and low density of p states at the Fermi level in AIB2 results in the 
dominance of the Fermi-contact mechanism. With the values of 0.7 (0.5) st./eV for the bare total (boron partial) 
DOS at the Fermi level for MgB2 and 1.7 eV for the effective Stoner exchange parameter, reasonable agreement is 
obtained with the experimental NSLR rates. In addition, the NSLR rate strongly depends on the concentration in 
the Mgi_a;Ala;B2 system, and the experimental studies of this dependence are highly desirable for the understanding 
of the anisotropy of boron p states. We note a relatively weak radial convergence of the NSLR rate in this system 
compared to the well studied d metals, as well as a common uncertainty related to the estimate of the enhancement 
factor. However, these uncertainties do not affect our qualitative conclusions. 

Similar results for NSLR in MgB2 were obtained in Ref. 139], where Ti on ^^Mg and the Knight shifts were also 
computed. We believe that the discrepancy in the Fermi-contact term is mainly due to the larger boron sphere 
radius used in Ref. [03. 
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Figure 9: Different contributions to the total ^^B NSLR in Mgi_2,AUB2 
thin sohd. Fermi-contact. 



Thick sohd hne, orbital; dashed, dipolar; 



4 Electron-phonon coupling 

Already in the first publication ||lg following the discovery of SC in MgB2 the strength of the EPC was estimated 
and a qualitative suggestion that MgB2 is a standard BCS superconductor was made. Measurements of the B- 
isotope effect |Q on Tc, tunneling |^, transport |Q, thermodynamic properties [Q, and the phonon density of 
states p4] confirm that MgB2 is most likely an electron-phonon mediated s-wave superconductor with intermediate 
or strong coupling. 

For the qualitative understanding below we use the simplest estimation (so-called rigid muffin-tin approxi- 
mation) with the formalism of Ref. Eq]. In this formalism the EPC constant A is proportional to the Hopfield 
parameter ry [Esl. To calculate r] one has to obtain the same complex off-diagonal density matrix at the Fermi 
level which enters the expressions for Ti (see Section 3). In Fig. nO we show the behavior of this parameter as 
a function of x in Mgi_a;Ala;B2 in the rigid band approximation. Changes in 77 are related to the changes in the 
electronic N , and even neglecting the changes in the phonon frequencies the general behavior is consistent with the 
experimentally observed trend for Tc in this alloy. The F-point frequencies of the E2g phonon mode obtained in our 
FLMTO calculation are 491 cm~^ in MgB2 and 956 cm~^ in AIB2. These frequencies are in reasonable agreement 
with other published data, as can be seen from Table 0. The biggest disagreement exists for the most sensitive and 
important E2g mode where different calculations produced results ranging from 470 to 665 cm~^. The hardening 
of the E2g mode with band filling should also contribute to the suppression of Tc with Al doping. Still, in spite of 
the qualitative agreement with experiment and low anisotropy of in-plane and out-of-plane contributions to 77, the 
theoretical Tc calculated in the rigid muffin-tin approximation is too low. 

It was suggested [g6| that holes at the top of the B-B bonding a bands may have strong coupling with the £^2g 
optical B-B bond-stretching mode. In Figs. |3p andHd it is shown how these 2D cylindrical sheets of the Fermi 
surface breathe with such distortion (see also Ref. ||25| ). The displacements in Figs. |^ and ^ are (c) O.OTas and 
(d) —O.OlaB (positive displacement corresponds to boron atoms moving towards the centers of hexagons). The 
internal cylinder (green in Fig. ^a) falls below Ep and is not seen. It is clearly seen that the bonding n band also 
experiences strong interaction with the E2g mode. 

State-of-the-art calculation of linear response |Q produced A — 0.87, tuiog — 504 cm^^ and ji* ~ 0.14. The 
corresponding Tc (from the solution of the Eliashberg equation on the real axis) is close to 40 K. The dominance of 
the a-a coupling via the optical bond-stretching mode was confirmed with the contributions from a and n bands 
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Table 3: Calculated phonon frequencies in MgB2 at F-point (in cm ^) 



Authors 


E2g 


Big 


A2u 


Ei„ 


Kortus M 




470 


690 


390 


320 


Yildirim [3 




486 


702 


402 


328 


Kunc 


46 


1 


535 


695 


400 


333 


Present 


491 


693 


389 


326 


Satta 


20 




665 


679 


419 


328 


Kong 


47 




585 


692 


401 


335 


Bohnen |4g 


536 


692 


394 


322 




Figure 10: Rigid band results for the Hopfield parameter 77 in Mgi_2;Al2;B2 obtained in non-spherical rigid-ASA 
approximation. Rigid band shift starting from MgB2 (solid line), and from AIB2 (dashed line). 
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to A amounting to 0.62 and 0.25, respectively, and the estimated uncertainty of about 6%. The final conclusion is 
that the unusually high Tj, is due to the large A caused by the presence of holes in the B-B bonding cr bands and 
by a relative softness of the E2g mode. It was also concluded that MgB2 is a clear case of an intermediate-coupling 



s-wave BCS superconductor. Other calculations pS, M, 48 confirmed these conclusions. In Ref. M strong non- 



linearity of EPC was found and compared with the measured neutron data. To illustrate the difference between 
MgB2 and AIB2, in Fig. O we show the phonon spectra from Ref. Eq] where A = 0.73 and 0.43 were found for 
these systems. 

In spite of the general agreement between the linear response and frozen phonon methods, it is unclear whether 
these perturbative techniques are adequate for such strong EPC. Unfortunately the detailed analysis of more 
sophisticated schemes of pairing WOl 10| lacks high accuracy of the linear response method even for this relatively 
simple material. 

Experimental estimation of EPC is not always consistent with the picture above. If we assume that 'bare' LDA 
N is essentially correct one can use several experiments to extract the strength of EPC. For instance, comparing 
7 from the specific heat measurements pO| with LDA 'bare' N according to the relation 7 = (l/3)7r^fc^7V(l + A) 
one would obtain A in the range of 0.53-0.62 (the LDA N is [18-20, 26] 0.70 ± 0.02 (eV-f.u.)-^). With such A 
an unusually low parameter /i* — 0.03 should be used to reproduce the experimentally observed T^. In addition. 



recent experiments 51 revealed no changes in the frequency of the i?2g phonon when the system goes through 
the SC transition. The analysis of the parameters {He (0)) /^T^ and AC (Tc) /"/Tc suggests [Q extremely weak 
coupling, and the thermoelectric power data |p3| cannot be fit with the LDA N. These controversies leave a 
question whether the electron pairing in MgB2 is phonon-mediated relatively open. We should stress that there 
is a large disagreement between different measurements. For example, the experimental data E3] for the specific 
heat 7 range from 1.1 to 5.5 mJ mol~^ K~^. Analyzing TV extracted from NMR data (see Section 3) one can 
find that the LDA partial boron N should be lowered by 10-15% to comply with the Ti measurements. With 
such correction for the total N the corrected A can easily be equal to 0.8. Thus, a lower N is required to make a 
satisfactory interpretation of these two experiments. At present the reasons for lowering of N are unclear. In any 
case the application of many-body techniques (GW or others) to this material should be valuable. 

5 Dielectric function and reflectivity 

Since optical properties reflect the electronic structure of a system, it is of interest to interpret the main features 
of the optical spectra in terms of the electronic transitions in the studied materials. Here we are dealing with 
the optical properties as expressed by the frequency-dependent dielectric function associated with direct interband 
transitions. 

The imaginary part e2(w) in the random phase approximation (RPA) with the LMTO basis set is p3| : 

4M = ^VE / dk\{kX\-^h^,\kX')\\f^il-f^,)Sie^,-el-hu;), 

^^ ^ ^ X X> JBZ 

where e\ and e\i are the eigenvalues for bands A and A', respectively, and f^ is the occupation number at zero 
temperature. For AIB2 structure, there are two independent components for E ±c and E [|c. 
The complex index of refraction is given by 

n^n + ik = e^l"^ = (ei + 162)^^^ , 

and the normal incidence reflectivity if defined by the Fresnel equation 

(n-l)2 + fc2 



R = 



{n + 1)2 + fc2 



We calculated the optical response for both polarizations. The tetrahedron method was used for the Brillouin 
zone integrations using 512 irreducible k points. The real part ei (w) was obtained from 62 (w) by the Kramers- 
Kronig relation. The e (w) functions were calculated up to 25 eV, which is well above the main peak of £2 {^)- Also, 
both in MgB2 and AIB2 there are partially filled bands, and therefore, one should take into account the intraband, 
or Drude, term. The plasma frequency in the Drude term is calculated from the average square of the Fermi velocity 
(see Section 2), and the scattering rate was obtained from the experimental data on resistivity in the normal state. 
This procedure is quite standard and was used for different metals including high-Tc superconductors p4, p5| . Since 
we have no accurate ellipsometric measurements of the dielectric function for these materials, our results should 
be considered as a prediction. 
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Figure 11: Theoretical phonon dispersion curves along the high-symmetry directions and phonon DOS for MgB2 
(top) and AIB2 (bottom). For details see Ref. Gsl. 
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Figure 12: Calculated e{Lo) for (a) MgB2 and (b) AIB2 with its band-to-band decomposition. 



Fig. ^ shows the interband contributions to £2(1^-') for i:^; < 11 eV for both materials and both polarizations. 
Also, we show the partial band-to-band contributions to the imaginary part of the dielectric tensor for all significant 
transitions. 

Both E ±c and E ||c components are very anisotropic and quite different for both materials. In MgB2, the 
function £2 (w) exhibits a structure with two broad peaks in E _Lc. The maximum at 2.5 eV corresponds to the 
transitions between the 4-th and 5-th bands mostly around the F-M direction at the Brillouin zone (see Fig. ||). 
The second maximum at about 8 eV is related to the 4<->6 transitions related to the H-A direction where these 
bands contain extended nearly parallel parts. For the polarization E ||c there is only one broad hump with several 
fine features which can be easily identified using the band-to-band transition analysis (Fig. |lja). The sharp peak 
at 5 eV is related to the 4^^5 transitions from narrow region around H-A direction where these bands are also 
nearly parallel. 

In AIB2 (Fig. |ljb) the two main peaks for the ab polarization are related to the 5<-^6 (3.8 eV peak) and 4^^6 
(7 eV peak) transitions, respectively. For the c polarization there is (like in MgB2) a single broad hump with the 
maximum at about 4.2 eV and with several noticeable features that are easily identifiable. 

There are two main differences between the MgB2 and AIB2 interband £2 (w) spectra. First, we notice that 
there is a low-energy peak (0.7 eV) in MgB2 which is absent in AIB2. This peak is related to the transitions 
between bands 3 and 4 in the H-A direction. The difference of energies for these two bands is very low in this 
direction. In MgB2 these bands cross the Fermi level while in AIB2 they are both below the Fermi energy, i.e., this 
peak does not appear. 

Second, the absolute values of £2 [to) are higher in AIB2 compared to MgB2 in all the considered energy interval 
and for both polarizations. The analysis of the joint density of states (which is equal to £2(^1;) when all matrix 
elements are set to unity) shows that it is related to larger matrix elements in AIB2. The reason for such behavior 
is unclear, and it is necessary to perform a detailed analysis of the wave functions for the bands surrounding the 
Fermi energy. 

The refiectivity spectra for MgB2 and AIB2 are shown in Fig. |l^. These spectra are anisotropic. For MgB2 
the anisotropy of the plasma frequency is not very high (Wp''= 6.5 eV, Lj'p=7.02 eV). However, the anisotropy of 
the reflectivity spectra is high even in the low-energy region due to interband transitions. For example, for the ab 
polarization there is a feature at ~ 1.3 eV corresponding to the peak of the interband part of £2(w) at ~ 0.7 eV 
(Fig. |3|). We believe that this structure can be observed in the refiectivity measurements. AIB2 does not have 
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Figure 13: Calculated reflectivity for MgB2 (left) and AIB2 
polarization. 



(right). Solid lines, c polarization; dotted lines, ab 



this low-energy peak, and therefore the behavior of the reflectivity at the energies below 2 eV is rather smooth. 

The studies of the general dielectric function e(k, w) and analysis of collective charge excitations for MgB2 
appeared recently in Ref. Q . 

6 Conclusion 

We showed that the electronic structure of Mg1_2.Ala.B2 has both similarity with and notable differences to GIC's. 
Its peculiar features include holes at the top of the a bands, coexistence of metallic and covalent bonding, and 
smooth disappearance of the 2D character in the electronic structure with Al doping. We emphasize how new NMR 
and optical experiments may provide a direct test of validity of these theoretical predictions. The presented results 
provide a basis both for further studies of normal and SC states of MgB2 and for the search of new SC compounds. 
Such studies are likely to be very complicated due to the fact that MgB2 seems to be a unique superconducting 
compound among the family of structurally similar materials with no or very low- Tc SC. It is unlikely that an 
analog to the high-Tc family of materials can be found in the case of MgB2. Nevertheless, the relative simplicity of 
this compound makes the development of new computational schemes with realistic electron-phonon coupling very 
attractive. 

The authors acknowledge useful discussions with F. Borsa, S. Bud'ko, P. Canfield, I. Mazin and N. Zein. This 
work was carried out at the Ames Laboratory, which is operated for the U.S. Department of Energy by Iowa State 
University under Contract No. W-7405-82. This work was supported by the Director for Energy Research, Offlce 
of Basic Energy Sciences of the U.S. Department of Energy. 
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